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Abstract 

The analytical representation of the longitudinal hadronic sho- 
wer development from the face of a calorimeter is presented and 
compared with experimental data. The suggested formula is partic- 
ularly useful at designing, testing and calibration of huge calorime- 
ter complex like in ATLAS at LHC. 



One of the important questions of hadron calorimetry is the question 
of the longitudinal development of hadronic showers [|IJ . For analysing of 
energy depositions in various depths of calorimeters it is necessary to have 
the analytical representation of longitudinal shower development from the 
face of calorimeter. Such formula is particularly useful at designing, testing 
and calibration of huge calorimeter complex like in ATLAS at LHC. 

There is well known parameterisation of longitudinal shower develop- 
ment from shower origin suggested by Bock et al. ||: 
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where X Q is the radiation length, A/ is the interaction length, a, b, d, w 
are parameters, k is the normalisation factor. In this work we use the 
values of parameters obtained in ||: a = 0.786 + 0.427 ■ InE, b = 0.29, 
d = 0.978, u = 1.032 - 0.159 • InE. 

In practice, the calorimeter longitudinal segmentation is unsumcient for 
precise determination of shower vertex position. Because of it is necessary 
to integrate the longitudinal profile from shower origin over the shower 
position: 
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where Xfj IS el coordinate of the shower vertex. 

In given work the analytical representation of the hadronic shower lon- 
gitudinal development from the calorimeter face has been derived. 

It turns out that 
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where iFi(a, 0, z) is the confluent hypergeometric function || and is 
the normalisation factor. The confluent hypergeometric function in the 
first term of formula (|3]) can be calculated by using the following relation: 

1 F 1 (l,a+ l,z) = az~ a e z -y(a,z), (4) 

where 7(0,2) is the incomplete gamma function ||. This function can 
be calculated by using the entry G106 in the program library ||. The 
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confluent hypergeometric function in the second term can be calculated by 
using the series development: 
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^(1, a + 1, -z) = 1 - + \ +.... (5) 

a + 1 (a + l)(a + 2) 

This series has the remarkable property: at a = it corresponds to the 
function e~ z . For real calorimeters with longitudinal size usually less then 
10 Xi the value of z — (1 — d)j^ does not exceed 0.2. Therefore it is 
sufficient only 3 terms in order to calculate the function iF\(l, a + 1, —z) 
with the precision of 0.1%. 
The normalisation factor is 

OO , . 

N = E beam / J (6) 
o 

where E beam is the beam energy. The integral in the denominator can be 
calculated by using the following relation [0: 

oo 

J z»e- Xz 1 F 1 (a 1 7 , kz) dz = T(v + l)\~ v ~ 1 2 F 1 (a, v + 1, 7 , ^) , (7) 
o 

where 2 Fi(a, /3,y,z) is the hypergeometric function ||, which in our case 
have the very simple form: 

2 F 1 (l,a+l,a + l,z) = j^—. (8) 

As a result we obtain 

AT = E beam 

A/ T(a) (w X b~ a + (l-w) \j d~ a ) ' v ' 

In Fig. m the calculations by formula (^) are shown and compared with 
the experimental data at 20 - 140 GeV obtained by using the conventional 
iron-scintillator calorimeter 01] and the one with innovative feature — the 
longitudinal orientation of the scintillating tiles |§. The good agreement 
is observed. 

So, now there is useful analytical formula which gives the possibil- 
ity immediately to obtain the longitudinal shower energy deposition from 
calorimeter face. 
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Figure 1: Longitudinal profiles of the hadron showers of 20 GeV (crosses), 
50 GeV (squares), 100 GeV (open circles) and 140 GeV (triangles) energies 
as a function of the longitudinal coordinate x in units Xj for conventional 
iron-scintillator calorimeter [|J and of 100 GeV (black circles) for tile iron- 
scintillator calorimeter ||. The solid lines are calculations by function (|3|) 
with parameters from J|. 
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